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CN ! I. INTRODUCTION 

^ , In a fundamental description of nucleon-nucleon (NN) interaction, the existence of the nucleon internal structure 
can not be ignored. The standard NN potentials are actually effective tools aiming to mimic a much more complicated 
interaction process, of which it is not even clear that it could be reduced to a potential problem. The NN system 
can be rigorously described only when starting from the underlying QCD theory for the nucleon constituents: quarks 
and gluons. This is however a very difficult task, which is just becoming accessible in lattice calculations [1-3], and 
that will be in any case limited for a long time to the two-nucleon system. In any attempt to describe the nuclear 
J> ' structure, one is thus obliged to rely on more or less phenomenological models. 

Since the nucleon size is comparable to the strong interaction range, the effects of its internal structure are expected 
to be considerable. In particular the NN interaction should be non-local, at least for small intcr-nucleon distances. 
In addition, we have no reason to believe that nuclear interaction is additive as the Coulomb one: the interaction 
between two nucleons may not be independent of the presence of a third one in their vicinity. Finally the interplay of 
, nucleon confinement and relatively large kinetic energies can generate - e.g. via virtual nucleon excitations ~ a rather 
■ strong energy-dependence in the interaction. Despite numerous studies devoted to this subject, we still do not have 
a clear understanding of the relative importance of these effects in NN force, specially concerning their influence on 
experimentally measurable quantities. 

This work investigates the consequences of using non-local NN forces in describing the A=3 and A=4 nuclear 
^ , systems. The locality of NN force, assumed in some of the so called realistic models [4, 5], is due more to numerical 
convenience than to convincing physical arguments. The two-nucleon experimental data, since they contain only 
on-shell physics, are successfully reproduced without including any energy dependence or non-locality in the NN 
^ ' force. However, they all suffer from the underbinding problem, i.e. two-nucleon interaction alone fails to reproduce 
k>( the nuclear binding energies, starting already from the simplest A=3 nuclei. Figure 1 shows the relative differences 
. between experimental and theoretical binding energies for He isotopes obtained with AV18 potential [4, 6, 7]. These 
d ' differences increase with the mass number A and vary from ~ 0.7 MeV in ■^He to ^ 10 MeV in the case of ^°He. 

The inclusion of non-local terms - like in Nijm 93 and Nijm I potentials [5] - does not remove this discrepancy 
[8, 9]. If in some cases, like in CD-Bonn [10] or in chiral models [11-13], they considerably improve 3- and 4-N binding 
energies, the improvement is still not sufficient to reproduce the experimental values. This underbinding is rather 
easily removed by means of three- nucleon forces (3NF). The existence of such forces is doubtless, but their strength 
depends on the NN partner in use and is determined only by fitting requirements. 

However the use of 3NF, to some extent, can be just a matter of taste. It has been shown in [14, 15] that two 
different, but phase-equivalent, two-body interactions arc related by a unitary, non-local, transformation. One thus 
could expect that a substantial part of 3- and multi-nuclcon forces could also be absorbed by non local terms. A 
considerable simplification would result if the bulk of experimental data could be described by only using two-body 
non-local interaction. In fact, the unique aim of any phenomenological model is to provide a satisfactory description 
of the experimental observables but it is worth reaching this aim by using the simplest possible approach. 

There exist already calculations reproducing the triton binding energy without making explicit use of 3NF. The 
first one was obtained by Gross and Stadler [16] using a relativistic equation and benefiting from some additional 
freedom in the off-shell vertex form factors. The complexity in using relativistic equations makes however difficult 
their extension to larger nuclear systems, or even to 3N scattering. 

A very promising result, which takes profit from non-locality in non relativistic nuclear models, has been obtained 
by Doleschall and collaborators [17-23]. In this series of papers, purely phenomenological non-local NN forces have 
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FIG. 1: "(Color online)" Comparison between the experimental and theoretical calculations with Argonne vl8 interaction in 
He isotopes. Results are taken from [6, 7]. 

been constructed, which were able to overcome the lack of binding energy in three-nucleon systems, namely "^H and 
■^He, without explicitly using 3NF and still reproducing 2N observables. The striking success of these models is closely 
related to the presence of small deuteron D-state probability, comparatively to local NN interaction. A difference 
which does not contradict the phenomenology and which follows from using two equivalent representations of the one 
and the same physical object [24]. In fact, non-locality of the Doleschall potential softens the short range repulsion 
of local NN models and can simulate part of effects due to the 6-quark structures as well as quark exchange between 
nucleons inside the nucleus. Local realistic interaction models, maybe with exception of Moscow [25] potential, 
artificially prohibits such effects by imposing a very strong short-range repulsion. 

Our work is an extension to the A=4 nuclei of the Doleschall pioneering studies. In particular, we would like to 
check whether such non-local interaction models remain successful or not when applied to the more complicated 4N 
systems. We will first provide results for ^He ground state and then extend the calculations to the n-'^H scattering. 
This system possesses a resonance at Ecm ~ 3 MeV, exhibiting different dynamical properties from those of bound 
states and testifying a failure of the conventional NN+3NF interaction models [26]. 

II. THEORETICAL TREATMENT 
A. FY Equations 

We describe the 4N system by using Faddeev-Yakubovski (FY) equations in configuration space [27-29]. Even 
though the major goal of FY formalism is a mathematically rigorous description of the continuum states, it turns out 
as well to be advantageous when dealing with bound state problem. The advantage lies in the natural decomposition of 
the wave function in terms of the so called Faddeev-Yakubovski components (FYC) which take benefit of the systems 
symmetry properties. These amplitudes have simpler structure than the wave function itself and are therefore easier 
to handle numerically. Four-particle systems require the use of two types of FYC, namely K and H . Asymptotes of 
components K incorporate 3-1-1 (see Fig. 2) channels, while components H contain asymptotes of 2-)-2 ones (see Fig. 2). 
By permuting the particles one can construct twelve different components of the type K and six components of the 
type H . The total wave function is simply a sum of these 18 FYC. 

In theoretical treatment of the problem we use the isospin formalism, i.e. we consider protons and neutrons as 
being degenerate states of the same particle - nucleon. For a system of identical particles, the FYC are not completely 
independent, being related by several straightforward symmetry relations. All the 18 FYC can be obtained by the 
action of the permutation operators on two of them, arbitrarily chosen, provided one is of type K and the other 
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FIG. 2: "(Color online)" Figure [a]: FY components K and H. Asymptotically, aa z ~* oo, components K describe 3+1 particle 
channels, whereas components H contain asymptotic states of 2+2 channels. Figure [b] illustrates the j-j coupling schemes used 
when developing FYC K and H into partial wave basis. 



of type H . We have chosen K = ^ffj 3 and H = '^\2- The four-body problem is solved by determining these two 
components, which satisfy the system of differential equations [29]: 

{E-Hq-V)K = V{P+ +P-)[{1 + Q)K + H] 

{E-Ho-V)H ^ Vi2p[{l + Q)K + H]. (1) 
P+, , P and Q are the permutation operators: 

Q = £^34 

P — P13P24 = P24:Pl3- 



Employing operators defined above, systems wave function is given by 

* = [1 + (1 + P+ + P-)Q] (1 + P+ + P-)K + (1 + P+ + P-){1 + P)H. (2) 



Components K and H are functions in configuration space, and depend also on the internal degrees of freedom of the 
individual particles (spins and isospins). The configuration space is provided by the position of the different particles, 
which we describe by using reduced relative coordinates. These coordinates differ from Jacobi coordinates usually 
employed in Classical Mechanics by factors depending on the particle masses. Use of such coordinates has several big 
advantages: first center of mass motion can be easily separated, then transition between two bases is equivalent to 
orthogonal transformation in space; finally, kinetic energy operator in this basis reduces to multidimensional 

Laplace operator in corresponding subspaces. Two principally different sets of reduced relative coordinates can be 
defined. One is associated with the components K = 5"' ■ ^: 



— -> _ /r) (m,+mj)mfc ^— > _ m^^rj+rnj^-. . . 



where mi and r g are respectively the mass and the position of the £-th particle. The coordinate set associated with 
the components H = ^fj, is defined by: 



2 . 



m -J^^iri-^k) (4) 



^. . > _ In ('"'+'"j)('"fc+'"i) /■ mfc'r'fc+mi r*i _ ™i r j+rrij r j x 

The functions K and H are expanded in the basis of partial angular momentum, spin and isospin variables, according 



to: 
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^ XiViZi 

Here Y°'{xi^yi, Zi) generalize tripolar harmonics containing spin, isospin and angular momentum variables. Func- 
tions J-'^{xi,yi, Zi) are so called partial amplitudes, being continuous in radial variables x,y and z. The label a 
represents the set of intermediate quantum numbers defined in coupling scheme, it includes as well the specification 
for the type of FY's component {K or H). We have used j — j couplings, represented in Fig. 2 [b], and expressed by: 




for components of K-type, and 




for the _ff-type components. Here Sj and U are the spin and isospin quantum numbers of the individual particles and 
{J'^ ,T) are, respectively, the total angular momentum, parity and isospin of the four-body system. Each amplitude 
J^"{xi,yi, Zi) is thus labelled by the a set of 12 quantum numbers a. The symmetry properties of the wave function 
in respect to exchange of two particles impose additional constraints. One should have (_y^+'^»=+*'" = £ for the 
amplitudes derived from any type of components {K or H), while for i?— type amplitudes additional constraint 
^_y!/+<^!/+tj/ = £ ig valid as well. Since we deal with nucleons (i.e. fermions) Pauli factor e is equal -1. The total 
parity tt is given by (^_y»=+'!'+'»^ ^ independently of the coupling scheme in use. 

By projecting each of the Eqs. (1) on its natural configuration space basis one obtains a system of coupled 
integrodifferential equations. In general one has an infinite number of coupled equations. Note that, contrary to the 
3N problem, the number of partial FY amplitudes is infinite even when the pair interaction is restricted to a finite 
number of partial waves. This divergence comes from the existence of additional degree of freedom 1^ in the expansion 
of the ii'-type components. Therefore we are obliged to make additional truncations in numerical calculations by 
taking into account only the most relevant amplitudes. 

B. Boundary conditions 

Equations (1) are not complete and should be complemented with the appropriate boundary conditions. Boundary 
conditions can be written in the Dirichlet form. First FY amplitudes, for bound as well as for scattering states, satisfy 
the regularity conditions: 

:r«(0, yi, Zi) = J^^ixi, 0, Zi) = r?{xi,yi, 0) = (8) 

For the bound state problem, the wave function is exponentially decreasing and therefore the regularity conditions 
can be completed by forcing the amplitudes to vanish at the borders of the hypercube [0, Xmax] x [0, i^ax] x 
[0,Zmax], i.e.: 

^ i (-^max5 2/i5 ^i) ~ ^i (^^ij ^^axj ^i) ~ ^i -^max) — (9) 

For the elastic scattering problem the boundary conditions are implemented by imposing at large values of z the 
asymptotic behavior of the solution. In case of N+NNN clastic scattering we impose at .^max the solution of the 3N 
problem for all the quantum numbers, corresponding to the open channel aa- 

° (xi, Z,„ax) = (a^i, Vi) (10) 

Functions f^^ixi^yi) are the Faddeev amplitudes obtained after solving corresponding 3N bound statc^ probkun. 
Indeed, below the first inelastic threshold, at large values of z, the solution of (1) factorizes into a bound st at solution 
of 3N Faddeev equations and a plane wave propagating in z direction with the momentum fc^^ = ^ ^ (-Ecm — E^n). 
One has: 
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There are two different ways to obtain the scattering observables. The easier one is to extract the scattering phases 
from the tail of the solution, namely taking logarithmic derivative of the open channel's K amplitude in the 
asymptotic region: 



tan(5= ' ' (H) 

This result can be independently verified by using integral representation of the phase shifts 



sinJ = -^j ^iT\^i'yiM''<^aZ){Vi4 + V24 + Vs4)^ {xi,yi,Zi)dV. (12) 

Here ^aa^\xi,yi) is a 3N bound state wave function composed by particles (1,2,3). This wave function is considered 
to be normalized to unity. Asymptotes of the wave function ^{xi, yt, Zi) is considered to have the same normalization 

as ^aa^\xi,yi), i.e. it tends to: 

^{xi, m, Zi) = ^'^^^^^ {xi, yi) {ka^Zi) + tan{S)hi^ {kaaZi)] ■ (13) 
A detailed discussion on these technical aspects can be found in [26, 30] . 



C. Numerical solution 



In order to solve the set of integro-differential equations - obtained when projecting eq.(l) in conjunction with the 
appropriate boundary conditions into partial wave basis - components J^" are expanded in terms of piecewise Hermite 
spline basis. 

J^nx, y.z)=Y. c':ikiSAx)Sk{y)Si{z) 

In this way, intcgro-diffcrcntial equations are converted into an equivalent linear algebra problem with unknown 
spline expansion coefficients c^^j to determine. In case of bound state problem eigenvalue-eigenvector problem is 
obtained: 



Ac = EBc, (14) 

where A and B are square matrices, while E and c are respectively unknown eigenvalue(s) and its eigenvector (s) to 
determine. In case of elastic scattering problem, a system of linear algebra equations is obtained: 



[A - E,mB] c = b (15) 

where b is an inhomogeneous term imposed by the boundary conditions eq.(lO). Numerical methods used for solving 
these large scale A'' ~ 10^ eigenvalue problems and linear systems are given in [26] . 



III. RESULTS AND DISCUSSION 



We have used in our calculations four different Doleschall potentials derived in references [17-23]. Hereafter, INOY96 
denotes the SB+SDA version of the potential defined in [17, 18]. It consists in short-range non local potentials in 

"'^So and ^SDi partial waves continued with a local Yukawa tail outside R=4 fm. The other partial waves arc taken 
from AV18. INOY03 denotes the IS version considered in Ref. [22]. It is an updated version of INOY96 which has 
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a smaller non locality range (R=2 fm) and provides a more accurate description of 2N observables. INOY04 and 
INOY04' are the two most recent versions [23] having the same ^Sq and ^SDi potentials as INOY03, completed with 
the newly defined non local potentials in P- and D- waves. As in the preceding models, higher partial waves are also 
taken from AV18. 

All the results have been obtained considering equal masses for neutrons and protons (m„ = nip = m) with 

^ = 41.47 MeV fm^. As mentioned in section II A we have used isospin formalism, furthermore assuming the total 
isospin quantum number T to be conserved. 

A. 3N system 

We will start with the presentation of onr results concerning 3N systems. Binding energies for '^H and ■^He nuclei 
are summarized in Table I. In order to control our accuracy, we have included in this table the results of AV18 [4] 
with and without Urbana IX three- nucleon force [31]. One can see that we are in close agreement with the benchmark 
calculations of Rcf. [32]. Our results concerning non-local potentials arc slightly different from those given in [22, 23]. 
The small deviations (« 15 keV) come from isospin breaking effects which were fully included in Doleschall calculations 
[33] while, as discussed in the above section, they were only approximately taken into account in ours. These isospin 
effects for AV18+UIX have also been evaluated in [32] and were found to be of about 5 kcV, i.e. three times less 
than in non-local models. In any case, the small differences related to isospin approximation can not overcast the 
main achievement of these non local interactions: the ability to reproduce experimental 3N binding energies without 
three-nuclcon force. One can also remark from Table I that, while the binding energies obtained with AV18-I-UIX 
are in good agreement with the experimental data, the value of AB = B^h — B^He is better reproduced by non local 
models. 



TABLE I: and ^He binding energies (in MeV) calculated with various versions of non-local Doleschall potentials and with 
AV18-I-UIX model. Results are compared to experimental values and previous calculations. AS denotes the difference between 
and *He binding energies 







^He 


AB 




this work 


other 


this work 


other 


this work 


other 


INOY96 


8.556 




7.882 




0.674 




INOY03 


8.497 


8.482 [22] 


7.734 


7.718 [22] 


0.763 


0.764 


INOY04 


8.476 




7.711 




0.765 




INOY04' 


8.464 


8.481 [23] 


7.704 


7.718 [23] 


0.760 


0.763 


AV18 


7.616 


7.618(2) [32] 


6.914 


6.917(2) [32] 


0.699 




AV18-FUIX 


8.473 


8.474(4) [32] 


7.739 


7.742(4) [32] 


0.734 




Exp. 


8.482 


7.718 


0.764 



The analysis of the 3N binding energies is shown in Table II. One can see that the major contribution to binding is 

due to the non-local short range interaction terms < Vni >■ Contributions from the local part - coming either from 
long range Yukawa tail or from higher partial waves - are marginal in INOY96 model and remain less than 15% in 
the other ones. 

Comparison between INOYs and AV18+UIX results is instructive. Both models provide similar energies and rms 
radii but they result from values of kinetic and potential energies which differ as much as ~ 50%. The introduction of 
a non-local interaction reduces the short range repulsion between nucleons and makes the potential well more shallow. 
Dciiteron is already obtained with the same binding energy and size as for AVI 8 model but its wave function does 
not have, at short distances, the sharp slope due to hard-core repulsion. On the other hand the weaker •^Si— ''Di 
coupling in non-local models generates a smaller contribution of D-state. All these effects reduce the average kinetic 
energy and favor stronger binding in the 3N system, which is more compact than deuteron. If the average size of 3N 
system between two models is only slightly different, the average kinetic energy of nucleons is sensibly smaller in case 
of non-local interactions. 

The relative contributions (algebraic values) of the most relevant Vnn partial waves to potential energy are given 
in Table III. The column labelled "others" denotes the contribution of all partial waves not listed in the table. One 
can remark the small role of P- waves. It was noticed long time ago that triton binding energy basically depends on 
^So and ^Si-^Di A/'A/'-interactions. Indeed, in absence of tensor force, the 3N ground state would have L=0 and 8=^ 
as conserved quantum numbers and in this case P-waves would not contribute at all. The '^Si-'^Di tensor coupling 
introduces an L = 2 admixture in the wave function. L=l state appears only in the second order and contributes less 
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TABLE II: Expectation values of kinetic (< T >) energies, rms radius {R = \/<r'^ >) and proton radius Vp corresponding to 
the binding energies of Tabic I. The values of the potential energy have been separated into contributions coming from the 
non-local < V„i > and local < Vj > terms of the potential. For AVIS+UIX model, < Vni > denotes the contribution of 3NF 
to potential energy. For ^He, the expectation values of Coulomb interaction have not been included. 





Model 


<T> (MeV) 


- <Vi> (MeV) 


-<V„i> (MeV) 


R (fm) 


rp (fm) 




INOY96 


34.24 


0.776 


42.02 


1.656 


1.561 




INOY03 


33.11 


5.551 


36.06 


1.664 


1.566 




INOY04 


33.01 


5.564 


35.92 


1.667 


1.567 




INOY04' 


32.97 


5.547 


35.89 


1.668 


1.568 




AV18 


46.71 


54.32 




1.770 


1.654 




AV18+IJIX 


51.28 


58.69 


1.140 


1.684 


1.584 




Exp. 










i.()() 


^He 


INOY96 


33.64 


0.777 


41.42 


1.684 


1.733 




INOY03 


32.33 


5.512 


35.21 


1.701 


1.752 




INOY04 


32.24 


5.510 


35.08 


1.703 


1.755 




INOY04' 


32.20 


5.525 


35.05 


1.704 


1.756 




AV18 


45.68 


53.30 




1.809 


1.867 




AV18+UIX 


50.22 


57.60 


1.095 


1.716 


1.767 




Exp. 










1.77 



than 0.1 %. NN P-wavcs start acting only in the second order as well, which explains their negligible contribution to 
3N binding energy, as shown in Table III. The reduction of the tensor force is also sizeable: ^Di-waves contribution 
are considerably smaller for Doleschall interactions than for AV18. 



TABLE III: Relative contributions of different Vjvjv partial waves to triton potential energy. 





3S1 


ISO 


3D1 


3P2 


3P0 


1D2 


3P1 


IPl 


3D2 


3D3 


Others 


INOY96 


57.94 


29.24 


12.78 


0.1658 


0.1506 


0.05664 


-0.3810 


-0.04528 


0.04939 


0.01505 


0.02518 


INOY03 


58.30 


30.31 


11.42 


0.1350 


0.1892 


0.04260 


-0.4214 


-0.05317 


0.06144 


0.01634 


0.02307 


INOY04 


58.35 


30.33 


11.38 


0.1409 


0.1579 


0.04843 


-0.4304 


-0.05495 


0.06303 


0.06460 


0.0113 


INOY04' 


58.40 


30.36 


11.37 


0.1478 


0.7854 


0.05061 


-0.4368 


-0.05711 


0.06527 


0.04333 


0.0089 


AV18 


45.00 


25.30 


28.97 


0.4466 


0.2272 


0.1840 


-0.3977 


-0.03244 


0.07952 


0.08769 


0.1290 



The calculated n-d scattering lengths are presented in Table IV. The quartet value (^a), corresponding to J'^ = |^ 
state, is independent of the interaction model in use, furthermore being in full agreement with the experimental one. 
This robustness is due to the strong Pauli repulsion, prohibiting two neutrons to get close; to each other. It follows 
that only ■'Si— '^Di V„p waves are important in describing J'^ = |^ state and still only through its, well controlled, 
long range part. Therefore, this state does not contain any off-shell physics and can be successfully described by any 
potential model, provided it reproduces the = 1+ np scattering observables. 

The integral representation of the phase shifts, eq. (12), is used to study the role of different Vjvjv partial waves. 
In Table V are given the relative contributions to this integral. Their sum, in algebraic values, is normalized to 100. 
Results for n-d doublet scattering length (^a) are presented in the upper half of the table and the quartet ones in the 
lower part. It can be seen that for J'^ = 1"^, NN waves other than ^Si contribute by less than 0.1%, which confirms 



TABLE IV: Neutron-deuteron (nd) scattering lengths (in fm) calculated using Doleschall potentials. 





^ttnd (fm) 


*a„d (fm) 


INOY96 


0.448 


6.34 


INOY03 


0.523 


6.34 


INOY04 


0.543 


6.:m 


INOYOi' 


0.553 




AV18 


1.26 


6.34 


AV18-FUIX 


0.595 


6.34 


Exp. 


0.65±0.04 


6.35±0.02 
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the statements above. The situation is different for ^a which results mainly from a cancellation between -'^Sq and ^Si 
and is more sensitive to higher NN partial waves, showing a deeper impact into the off-shell physics. Due to the 
smallness of ^a, all the interaction effects have to be taken into account very accurately. In particular, its value is 
very sensitive to the electromagnetic (e.m.) interaction terms. The differences between INOY predictions and the 
experiment can therefore be caused by the absence of e.m. terms in these models. On the contrary e.m. corrections 
were properly included in AV18+UIX results. In any case the small discrepancy with data has no consequences in 
phenomenology, since ^ a is by one order of magnitude smaller than "^a and its relative contribution to the scattering 
cross sections is negligible. 



TABLE V: Relative contributions of different NN interaction waves in n-d integral scattering lengths. Doublet value is in the 
upper half of the table and quartet in the lower 





3S1 


ISO 


3D1 


3P2 


3P0 


1D2 


3P1 


IPl 


3D2 


3D3 


Other. 


INOY96 


-685.3 


800.3 


-20.14 


2.664 


-1.240 


-5.097 


-11.97 


20.07 


-0.7484 


0.6967 


0.7998 


INOY03 


-572.6 


685.2 


-16.73 


1.955 


0.2399 


-4.523 


-10.95 


16.73 


-0.6048 


0.3972 


0.8616 


INOY04 


-549.6 


662.4 


-16.28 


1.997 


-0.3862 


-4.257 


-10.31 


15.66 


-0.5025 


0.3749 


0.8932 


INOY04' 


-538.5 


650.8 


-15.92 


2.127 


-1.076 


-4.167 


-9.308 


15.25 


-0.4694 


0.3563 


0.8825 


AV18 


-195.5 


293.5 


-4.283 


4.779 


0.6261 


-0.3768 


-5.944 


6.190 


-0.3316 


0.9574 


0.3281 


INOY96 


100.1 


-0.0297 


0.0183 


-0.2708 


-0.5036 


-0.1510 


0.8177 


0.0190 


-0.1644 


0.0076 


0.0312 


INOY03 


100.0 


-0.0288 


0.0199 


-0.2694 


-0.4640 


-0.1529 


0.8180 


0.0192 


-0.1641 


0.00742 


0.0315 


INOY04 


100.1 


-0.0283 


0.0197 


-0.2707 


-0.4938 


-0.1515 


0.8343 


0.0189 


-0.1653 


-0.0006 


0.0480 


INOY04' 


99.98 


-0.0279 


0.0195 


-0.2675 


-0.4973 


-0.1503 


0.9043 


0.0188 


-0.1651 


-0.0004 


0.0479 


AV18 


100.1 


-0.0400 


0.0097 


-0.2672 


-0.4870 


-0.2577 


0.8081 


0.0252 


-0.1691 


0.0096 


0.0361 



B. 4N system 

Our results concerning the a-particle binding energy are displayed in Table VI. Two series of calculations were 
performed, including (upper half of the table) and neglecting (lower half) Coulomb repulsion between protons. This 
latter interaction was provided by Argonne group in their AV18 code [4] and takes into account proton finite size 
effects. 



TABLE VL Binding energy B (in MeV) and rms radius R (in fm) for *He ground state obtained with Doleschall and AV18-I-UIX 

potentials. The lower part contains Coulomb force. Energies presented in the two last lines of the table respectively for AV18 
and AV18-I-UIX models have been taken from reference [8, 34], whereas rms radius from reference [6]. 



Pot. 


< T > 


-<V > 


B 


R 


INOY96 


72.80 


103.8 


31.00 


1.353 


INOY03 


69.89 


99.94 


30.04 


1.369 


INOY04 


69.49 


99.41 


29.91 


1.372 


INOY04' 


69.46 


99.36 


29.88 


1.372 


AV18 


98.69 


123.6 


24.95 


1.511 


Pot. 


< T > 


-<V > 


- <E> 


R 


INOY96 


72.45 


102.7 


30.19 


1.358 


INOY03 


69.54 


98.79 


29.24 


1.373 


INOY04 


69.14 


98.62 


29.11 


1.377 


INOY04' 


69.11 


98.19 


29.09 


1.376 


AV18 


97.77 


122.1 


24.22 


1.516 




97.80 


122.0 


24.23 [8, 34] 




AV18-FUIX 


113.2 


141.7 


28.50 [8, 34] 


1.44 [6] 


Exp 






28.30 


1.47 



Calculations have been done by considering isospin averaged pair interaction, i.e.: 

^tit2 = Pnn{tl-,t2)Vnn + -Ppp(il , ^2) V^p + -Pnp (* 1 ; * 2 ) V„p 

where (ii, 12) arc the isospin quantum numbers of FY amplitudes in equation (5). They respectively represent (ij;, T3) 
for if- type amplitudes and {tx,ty) for if -type. PnniPpp and P„p are the probabilities of finding respectively nn, pp 
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and np pairs in a given isospin state. Note that, since the number of protons and neutrons in a-particle is equal, one 

has P„„(ti,t2) = Ppp{ti,t2)- 

As in 3N calculations, wo have neglected isospin breaking effects, considering a-particle as a pure T = state. 
Contributions of T = 1,2 admixture were calculated for AV18 and AV18+UIX models in [8] and found to be as small 
as 10 keV. Results for 3N system presented in last section showed that Doleschall non-local models are more sensible 
to isospin breaking: they account for «15 keV in ''H compared to w5 keV in AV18+UIX [32]. In any case, for the 
alpha particle these effects should not exceed some 50 keV and will not affect the physics discussed below. Notice 
also that Coulomb corrections obtained by non-local models exceed by 70 keV those obtained by AVIS, due to the 
different rms radii they give. 

As mentioned in section (II A), FY calculations have been performed in the j — j coupling scheme. The following 

truncations in the partial wave expansion of amplitudes were used: (i) Ynn waves limited to Za; < 3 but always in- 
cluding tensor-coupled partners, i.e. involving the set ^^o,^ SDi, iPi,^ Pq,^ PF2? Pi, ^D^? DGa,^ D2, ^^3,^ FH^,^ F3 
and (ii) Ix + ly + lz < 10. 

Convergence was studied as a function of jyz=max(jy,jz) for K-like components and jyz=max{jy,lz) for H-like, 
starting with jy^ < 1. In Table VII we present the a-particle binding energy results for INOY04' and AV18 models 
respcictivcly. The convcirgcnicc! is rathcir smooth, cixcept when passing from jy^ < 5 to jy ^ < 6. We think this is an 
artefact of our truncation procedure which keeps fixed the basis set in the x-coordinate. Note that the agreement 
between our results for AV18 potential and those, given in reference [8, 34], is very good. From results displayed 
on Table VII, as well as from analogous convergence patterns seen in 3N calculations, we conclude that Doleschall 
potentials converge more rapidly than AV18. This is probably due to their weaker tensor force, resulting into wave 
functions with stronger spherical symmetry. 



TABLE VII: Results of a-particle binding energy (in MeV) for INOY04' and AV18 models with Vatjv interaction limited 
to Ix < 3 (sec text) and partial wave basis limited to Ix + ly + Iz < 10. The convergence was searched as a function of 
jyz—inayi{jy,jz) for K-like components and jyz=ma,x{jy,jz) for H-like. Last line, denoted by *, contains additional calculations 
with NN interaction waves up to jx < 6 and Ix + ly + Iz < 12. 



Jyz 



INOY04' AV18 



28.094 
28.661 
28.967 

28.971 23.897 
28.974 23.920 

29.084 24.233 

29.085 24.226 
29.085 24.223 



To our opinion the main conclusion of the results displayed in Table VI is the possibility offered by the INOY 
models to provide a satisfactory description of A=4 nuclei in terms of two-body forces alone, as it is already the 
case for A=:2 and 3. One can argue that this agreement is not yet fully realized in their present version, for they all 
slightly overbind the experimental value: the most favorable version (INOY04') still exceeds by 0.79 MeV the "'He 
binding energy. One should however remark that this result is obtained without adjusting any additional parameter 
with respect to A=3. On the other hand, the difference between INOY96 and INOY04' due essentially to different 
parameterizations of their non local short range parts - is 1.1 MeV. It seems thus possible, by a finer tuning, to reach 
an even more precise description of A=4 in a next generation of potentials. If they are not contradicted by other 
aspects of the phenomenology, INOY models offer an alternative description permitting to avoid three-nucleon forces. 

Results of Table VI have been gathered in a Tjon plot - see Figure 3 - which displays the correlation between and 

''He binding energies for various NN potentials. One can see that, due to the small overbinding of a-particle, INOY 
results (diamond symbols) are outside the line formed by realistic local model predictions and, except for INOY96, 
are almost superimposed to CD Bonn -|- TM value. 

INOY03, INOY04 and INOY04' models, which differ in their P-wave structure, give very close results while INOY96, 
which has a different non-local S-wave structure, falls out further apart. This indicates that the S-waves are the key 
point in binding a-particle and in order to improve the agreement with the experimental value a better tuning in the 
^Si-^Di and ^So could be helpful. 

Proton rms radii predicted by INOY potentials deserve some comments. One can sec already in 3N systems (see 
Tab. II) that they are slightly smaller than the experimental ones. For a-particle we have only calculated average 
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rms of nucleons, without making distinction between neutrons and protons. The real value of protons rms should be 
slightly larger. However, since Coulomb interaction has very small effect on the alpha particle wave function (rms 
radii calculated by taking Coulomb interaction into and without account differ only in fourth digit) this value can 
not differ by more than 0.5%. One can therefore see that protons rms provided by INOY model are already by 6% 
smaller than the experimental one, compared to 1-2% in 3N complex. This fact clearly demonstrates that INOY 
interactions are too soft, resulting into a faster condensation of the nuclear matter. In order to improve the agreement 
one should try to increase the short range repulsion between the nucleons. This would inevitably imply a reduction 
of 3N binding energies, although this reduction is not necessarily very large. If one chooses the ITF ■^Si-'^Di and 
ISA ^So potential versions from [21], which give the smallest probabilities to find nucleons close to each other in 2N 
systems, the resulting triton underbinding will be of only ~50 keV. 



T ' I ' I ' I ' I ' r 




H 1 1 1 1 1 , 1 1 1 1 1 1 

7.6 7.8 8.0 8.2 8.4 8.6 

B (^H) [MeV] 

FIG. 3: "(Color online)" Tjon-line for the local and non local NN potentials. 

Let us finally consider the n+^H elastic scattering. This is, in principle, the simplest 4N reaction, being almost a 
pure T — 1 isospin state and free of Coulomb interaction. However, its simplicity is only apparent. In fact, n+^H 
is a system with very large neutron excess (77 = = 0.5). The only stable nucleus having a neutron excess 

equally large is ^He. Let us remind that AV18+UIX hamiltonian, faces increasing difficulty when describing neutron 
rich nuclei. The more elaborate 3NF model, namely Illinois [6], even though being able to improve the agreement 
with experimental data, still suffers from similar discrepancies [26]. In addition, the n+'^H system contains several 
nearthreshold resonances of negative parity (see Fig 4) , which strongly affect the scattering observables near Ecm ~ 3 
MeV. Their internal dynamics is richer than for bound states and it is not clear that they can be described by using 
the same recipes than the one used for solving the underbinding problem. The description of the n+'^H cross sections 
in the resonance region is therefore a very challenging task for nucleon-nucleon interaction models. 

We have performed extensive calculations of the n+^H scattering states only by using the INOY04 potential. The 
model dependency of the results was checked at Ecm—0 and 3 MeV with INOY04'. We present in Table VIII the 
calculated singlet (ao+) and triplet (ai+) scattering lengths together with the deduced coherent value 

Oc = ^ (ao+ + 3ai+) (16) 

and the zero energy cross section 



(17) 



11 



> 

CD 



3760- 



3756- 



3752 



3748- 



nnnp 
dnn 



3lj 

n H 



r 

0' 

r 

2" 



FIG. 4: "(Color online)" Negative parity resonance above the n+t threshold 



TABLE VIII: 


n+^H singlet (J''=0+), triplet (J'' = l+) 


and coherent scattering 


lengths (in fm) along with zero 


energy cross 


sections. Different model results are compared with the 


experimental data. 






Pot. 


%)+ (fm) 


ai+ (fm) 


a.c (fm) 


a (fm") 


INOY04 


4.00 


3.52 


3.64 


166.5 


INOY04' 


4.00 


3.52 


3.64 


166.8 


AV18 


4.27 


3.71 


3.85 


187.0 


AV18+UIX 


4.04 


3.60 


3.71 


173.4 


Experimental 


3.70±0.62 


3.70±0.21 


3.82±0.Q7 [39] 


170±3[38] 




4.98±0.29 


3.13±0.11 


3.59±0.Q2 [40] 






2.10±0.31 


4.05±0.09 








4.45±0.10 


3.32±0.02 


3.607±0.017 [41] 





Results for AV18 and AV18+UIX models have also been obtained and agree at 1% level with those given in Ref. [35]. 

The J'^ = 0^ and 1"*" positive parity states, determining the low energy behavior of the n+^H cross section, do not 
have any jS-matrix singularity, except the triton bound state threshold. It is therefore not surprising that the n+'^H 
scattering lengths are found to be correlated with 3N binding energy, in a similar way as n + d doublet scattering 
length is [36]. This is the reason why realistic local interaction models, providing too low 3N binding energies, 
overestimate n+^H zero energy cross sections. Once triton binding energy is corrected, for instance by implementing 
3NF, a value close to the experimental one is automatically obtained. From Table VIII it can be seen that Doleschall 
potential agrees with the lower bound of experimentally measured zero energy cross section, whereas AV18+UIX 
model coincides with its upper bound. The zero energy scattering cross section is thus fairly well reproduced. 

The situation with scattering lengths looks more precarious. The values found in literature are hardly compatible 
with each other [37] as it can be seen in Table VIII. The usual way to get is to express them in terms of the 
measured quantities Oc and cr{0), by reversing relations (16) and (17). This procedure, represented in Fig. 5, is 
numerically unstable. Indeed, once cr(0) is fixed, the domain of permitted ai+ and ao+ values is given by the ellipse 
of eq. (17) in the (ao+,ai+) plane. Since there are uncertainties in (t(0), the permitted values of scattering lengths 
are trapped in-between two ellipsis (dot curves in Fig. 5). On the other hand, each measurement of Oc restricts ai+ 
and ao+ values to lie on a straight line which spreads into a band due to experimental errors (see Figure) . The lower 
band displayed in Figure 5 follows from the R-matrix analysis result Oc — 3.607 ± 0.017 fm [41], while the upper 
one comes from the experimental measurement Oc — 3.82 ± 0.07 fm from [39]. By assuming an exact value of Cc, 
e.g. Qc = 3.624 fm given by the top of the lower band, the present - though small - experimental error in (t(0) 
leads to two sets of solutions which spread over a wide range: (i) «□+ = [4.31 — 5.00], ai+ = [3.16 — 3.40] and (ii) 
ao+ — [2.25 — 2.94], ai+ — [3.85 — 4.08] fm. This example illustrates the difficulty of extracting reliable values of ao+ 
and ai+ . The accurate determination of would require to gain one order of magnitude in measuring both (t(0) and 
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2 3 4 5 

a,, (fm) 



FIG. 5: "(Color online)" Extraction procedure for n+^H singlet (ao+) and triplet (aj+) scattering lengths from measurements 
of zero energy cross section (elliptic band) [38] and coherent scattering length (linear bands) [39-41] . The values of ai axe given 
by the intersection of these two curves. Band widths are related to experimental errors and, even being small, they make their 
determination very unstable. 

As it can be seen also from Figure 5, the coherent scattering length value Uc ~ 3.82 ± 0.07 fm of reference [39] is 
in evident disagreement with the experimentally measured zero energy cross sections, since it does not intersect cr(0) 
ellipsis. In this respect, the more recent values Oc = 3.607 ± 0.017 fm [41] and Oc = 3.59 ± 0.02 fm [40] are more 
reliable. Doleschall non-local potential provides = 3.63 fm, one standard deviation from these measurements, and 
seems to be more compatible with data than AVIS+UIX model. Fig. 5 suggests also that the real value of the zero 
energy cross section should coincide with the lower bound of the experimental result. 

The success in describing n+'^H scattering lengths by Doleschall potential is visible at slightly higher energies as 
well. In Figure 6 we present our calculated elastic cross section for the scattering energies in the n+^H center of 
mass energy range from to 3 MeV. Doleschall potential reproduces experimental cross sections near its minima 
at Ecm « 0.4 MeV. In this region both MT I-III - the only potential known to us being capable to reproduce the 
resonant region [43] - and AV18+UIX overestimate the experimental value. 

In previous works [30, 37, 42, 44, 45] we pointed out that local realistic interaction models underestimate the 
cross sections near the resonance peak, Ecm —5 MeV. At that time, calculations had been however performed with 
a limited number of partial waves and the failure was attributed in Ref. [47] to a lack of convergence. Recently we 
have considerably increased our basis set and have shown that the disagreement is indeed a consequence of nuclear 
models [26, 46]. The convergence of our present results is shown in Table IX, following the same truncation criteria 
as for a-particle (see Table VII). The number of FY partial amplitudes involved in n+t scattering calculations is 
considerably larger than for a pure 0"^ bound state and we have not been able to go in the PWB as far than in Table 
VII. One can however remark that results displayed on Table IX converge pretty well, and provide at least three digit 
accuracy. 

Implementation of 3NF is just able to improve zero-energy cross sections and is not efficient at the resonance 
energies. Doleschall non-local potentials seem to suffer from a similar defect: the phase shifts obtained using INOY04 



13 



3.0 



2.5- 



2.0- 



1.5- 



1.0- 



0.5- 



0.0 




0.1 



--MT l-lll 

• AV.18+UIX 

♦ Doleschall 



E (MeV) 



FIG. 6: "(Color online)" Comparison between experimental and theoretical n-'^H total cross section calculated with several 
local and non-local NN potentials. 



TABLE IX: Convergence of n+t scattering lengths and selected phaseshifts at Ecm=2.625 MeV for INOY04 model. Corre- 
sponding mixing parameters are given in parentheses. 





a+ (fm) 


a+ (fm) 


5(1+) (dcg) 


S{1-) (deg) 


1 


3.889 


3.609 






2 


3.995 


3.508 


-56.13 -0.757 (0.803) 


21.32 39.54 (-42.05) 


3 


3.995 


3.513 


-56.12 -0.759 (0.803) 


21.39 39.74 (-43.06) 


4 


3.995 


3.515 


-56.12 -0.759 (0.803) 


21.39 39.76 (-43.16) 



are even slightly smaller in their absolute value (except for 2~ state) than those obtained with local potentials and the 
total cross section is slightly worse. In fact the reduction of positive-parity phase shifts is a consequence of improving 
triton binding energies. As was previously discussed, n-t-"^!! scattering lengths are linearly correlated with the triton 
binding energy. Whatever the way one uses to increase triton binding, by means of non-local interaction or 3NF, 
the final result will inevitably be a reduction of 0+ and 1+ n+^H scattering lengths and low energy phase shifts (in 
absolute value). It turns out that by the same way, we reduce in absolute value the 0~ and 1~ phase shifts. Only 2~ 
phases are slightly increased in both AV18-I-UIX and Doleschall models. Thus we have a real puzzle for the interaction 
models: on one hand they have to reduce low energy cross sections, while on the other hand cross sections in the 
resonance region should be significantly increased. The fact that all realistic interactions systematically suffer in the 
resonance region let us believe that the underlying reason of this disagreement is not related to the non-locality or to 
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3NF effects. The observed discrepancies have different background than the underbinding problem. Where does this 
failure comes from? 

When analyzing n+^H cross sections in the resonance region, one should first recall their origin. They are negative 
parity states and their symmetry is consequently different from the positive parity ones, which are dominated by 
S-waves. The good agreement of the scattering lengths provided by Doleschall potential as well as its success in 
reproducing the low energy cross section minima, makes us believe that the positive parity phases are quite well 
reproduced in the resonance region as well. On the other hand, negative parity phase shifts should be rather far from 
reality, causing a disagreement with the experimental data. 



TABLE X: Relative contributions of different NN interaction waves in n-^H integral scattering lengths (second half of the 
table) . 





r 


Ecm (MeV) 


ISO 


3S1 


IPl 


3P0 


3P1 


3P2 


1D2 


3D1 


3D2 


3D3 


others 


INOY04' 


0+ 


0.0 


75.95 


22.83 


3.588 


-1.301 


1.942 


-0.8112 


-0.8661 


-1.433 


0.1145 


0.006131 


0.0188 


INOY04' 


0+ 


3.0 


79.79 


19.68 


3.413 


-1.046 


1.913 


-0.5796 


-0.5796 


-1.837 


0.0745 


0.005113 


-0.0803 


INOY04 


0- 


3.0 


61.93 


67.82 


-0.3569 


32.97 


-26.58 


2.190 


1.897 


-40.41 


0.6489 


-0.6528 


0.5551 


INOY04' 


0- 


3.0 


64.16 


70.14 


-0.3769 


32.83 


-29.71 


2.305 


1.987 


-41.90 


0.6723 


-0.6919 


0.5750 


INOY04' 


2- 


3.0 


39.75 


54.90 


-0.1640 


1.063 


-6.893 


14.29 


0.0970 


-3.182 


0.1476 


.0001 


0.0092 



To understand the possible source of such a disagreement, we have calculated - as for nd case - the relative 

contribution of the different NN partial waves in the integral expression of the phase shifts (12). The obtained results 
are summarized in Table X. First two rows correspond to the 0+ at E=0 and E=3 MeV. One can see that positive 
parity states are completely controlled by the interaction in ^SDi and ^Sq waves, at zero energies as well as at Ecm=3 
MeV, close to resonance peak. The role of higher partial waves is marginal. The sitiiation changes dramatically in 
negative parity states (values on 3-5 rows). Contribution of P-wave interactions becomes comparable to S-wave. On 
the other hand, the non triviality of physics in the resonance region is reflected by the strong compensation of different 
P-waves as well as ^Di and components in ^SDi channel. In addition, different P-waves dominate in different 
states: in 2~ state, the ^P2 waves are the most relevant, whereas ^Pq is almost negligible; in 0~ state, the ^Pq wave 
has the largest contribution, while 3P2 fades away. 

Finally, we would like to comment that all the observables where NN P-waves are contributing, have tendency to 
disagree with the experimental data. A small disagreement can already be seen in n-d doublet scattering lengths 
(Table IV), whereas triton, described by the same quantum numbers but where P-waves arc negligible, is perfectly 
reproduced. Other examples could be the 3N analyzing powers [49, 50], as well as increasing discrepancy when 
describing binding energies of neutron rich nuclei (see Figure 1). One should also recall that P-waves in most of the 
NN interactions are tuned on n-p and p-p data. Moreover, p-p P-waves are overcasted by Coulomb repulsion, while 
n-n P-waves are not directly controlled by experiment at all. This study suggests that CSB and CIB effects can be 
sizeable in n-n P-waves and provide a possible explanation for the disagreement observed in n-^H resonance region. 

IV. CONCLUSIONS 

During the last decade, a series of non-local NN potentials has been developed by Doleschall and collaborators and 
were found to provide an overall satisfactory description of the 2N and 3N system. If they left unsolved some of the 
theoretical nuclear problems - like the so called Ay puzzle - they constitute an undoubted success for their ability to 
reproduce the experimental binding energy of and ^He nuclei without adding three-nucleon forces. 

In this work, we have examined the possibilities for these non-local potentials to describe the 4N system as well. 
This system is a cornerstone in the nuclear ab-initio calculations and a crucial test for the nuclear models. Not 
only because therein the underbinding problem manifests in its full strength, but also because of the rich variety of 

scattering states it possesses. 

We have found that non-local NN models, could well provide 3N and 4N binding energies in agreement with the 
experimental data without making explicit use of three-nucleon forces. They offer an alternative solution to cope with 
the nuclear underbinding problem, than the one offered by the local plus 3NF philosophy. 

In their present form, they overbind "^He by some 0.7 MeV, a discrepancy much smaller that all the existing models 
and which reverses the lack of binding observed in most realistic potentials. The fluctuations in the predictions 
between different versions of the non-local models suggest that a finer parametrization of their internal non-local part 
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could be enough to make them fully successful in that point. The current version seems to be too soft in the short 
range region, thus giving slightly too small rms radii of light nuclei. 

By calculating n+^H scattering states, we have found that Dolcschall models arc also very encouraging in describing 
the low energy parameters, providing an even better description than Nijm II or AV18+UIX. However they fail also 
in reproducing the elastic cross section few McV above, in the resonance region. In a series of preceding works 
[26, 42, 44, 46] we have shown that local realistic interactions, even implemented with 3NF, underestimate the cross 
sections at the resonance peak Ecm= 3 MeV. Unfortunately, the non local models do not solve this problem. On the 
contrary, they even provide slightly smaller values of the cross section. Wc believe than the reason of this failure is 
common to all realistic models and lies in the nucleon-nucleon P- waves themselves. The analysis of their contribution 
shows that, contrary to ^He binding energy, they play a crucial role in the n+^H cross sections. If the n-p P- waves 
seem to be well controlled by the experimental data, one still has a relative freedom in the n-n ones to improve the 
description. 
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